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An electronic excitation of a molecule by a sequence of two femtosecond phase-locked laser pulses is 
considered. In this case the interference between the vibrational wave packets induced by each of the 
subpulses within a single molecule takes place. It is shown that due to the dynamical squeezing effect 
of a molecular vibrational state the interference of the vibrational wave packets allows one to measure 
the duration of a femtosecond laser pulse. This can be achieved experimentally by measuring the 
dependence of the integral fluorescence of the excited molecule on the delay time between the 
subpulses. The interference can lead to a sharp peak (or to a down-fall) in that dependence, the width 
of which is equal to the duration of the laser pulse. It is shown that finite temperature of the medium 
is favorable for such an experiment. 


Recently a great interest has been shown to the study of spatially localized vibrational wave 
packets in molecules induced by ultrashort laser pulses. The time evolution of the mean position of 
such a wave packet corresponds to the nearly classical nuclei motion, and can be observed in the 
pump and probe spectroscopic optical experiments [1-4]. The idea of controlling the chemical 
reactions due to the possibility of the nearly classical nuclei motion has been widely discussed and 
experimentally verified [5-10]. 

An important characteristic feature of the vibrational wave packet is its spatial extent. In the case 
of a harmonic nuclear potential the spatial properties of the wave packet are closely connected 
with the phonon statistics and can be treated by the methods developed in quantum optics. In a 
series of previous publications [1 1-16] we have studied the statistical properties of the vibrational 
states appearing due to the excitation of Franck-Condon transitions by the transform-limited light 
pulses of finite duration. The problem was to learn how is it possible, by varying the characteristics 
of a laser pulse (i.e. its amplitude and phase modulation), to excite the molecular vibrations with 
the given quantum fluctuations of the conjugated variables. In particular, it was theoretically 
predicted [11-13] and experimentally verified [17] that, when the spectral width of the exciting pulse 
is smaller than that of the absorption band, there appears a molecule in a squeezed vibrational 
state with a reduced quantum uncertainty of the nuclei position. The mechanism of squeezing, 
which arises here, is of the dynamical nature and can be explained as the result of quantum 
interference in the phase space of the molecular vibrations [13]. An application of the dynamical 
squeezing effect to the problem of wave packets optimal shaping for the control of the chemical 
reactions was considered in [18-21]. 

Our goal now is to show how this dynamical squeezing effect can be put to use in another way 
for the duration measurement of a femtosecond laser pulse. Following [22], we shall consider the 
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excitation of a Franck-Condon transition in a molecule by a sequence of two femtosecond phase- 
locked Gaussian-shaped pulses: 

E(/) = ^y-{ ex P - y( / + T /l\ + « ex P jexp(-/Q/) + c.r., (1) 

separated from each other by the time interval T . Here a is a possible additional complex 
coefficient between the pulses: a = |a|exp(/^). As it was demonstrated experimentally in [22] it is 
possible to change the time interval T , keeping constant the value of phase (p . 

We shall assume that due to the electronic transition only a spatial shift of the harmonic nuclear 
potentials occurs in the molecule. The adiabatic Hamiltonians describing molecular vibrations in 

the initial (/) and excited (?) electronic states have the form 

* 2 2 

p mco * . 

H ‘ =€i+ 2^ + — (q + qi) ' (2) 

P 2 mo}1 -o 

"' =£ ' + ^ + ~»' <3) 

where e ie are the electronic energy levels. It is convenient to rewrite Eqs. (2) and (3) expressing q 

and p through the phonon creation and annihilation operators: 

< 4 > 

Then 

mcojq, tico / \ / \ \ma> 

H, = e, + — - — + ~{b b + bb^) + ha^b' +b), g= q, y— , (5) 

H e = £ e +—^~(b + b + bb + ). (6) 

A non-dimensional coupling constant g , which appears in Eq. (4), is equal to the ratio of the 
Franck-Condon shift to the amplitude of zero vibrations. The Hamiltonian (5) can be diagonalized 
by means of the unitary displacement operator 

D = exp[-g(Zf - />)], (7) 

D * H t D = €j + ~T~{b b + bb + ). ( 8 ) 

It is easily seen from Eq. (8) that the ground state vibrational wave function of the molecule in 
the initial electronic state has the form 

|O f >= Z3|0>, (9) 

where |0) is the phonon vacuum, which coincides with the ground vibration state of the 
electronically excited molecule. Let us note that the displacement operator acting on the vacuum 
state gives the coherent state. Then Eq. (9) simply means that the ground state wave function of the 
harmonic oscillator being placed in the shifted potential is a coherent state. 

The case of zero temperature. We shall assume that initially at / = -ao the molecule was prepared 
in the state |/)|d> ( ), where |/) is the electronic wave function of the unexcited molecule. After the 
resonance interaction with the field ( 1 ) the wave function of the molecule will have the form 


|>P,r>= If.^l «>,.»>-' 




<< 1. 
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Eq. (10) is derived in the first-order perturbation theory. Here d (> , is the dipole matrix element of 
the electronic transition, which is assumed to be independent of the vibrational coordinate in 
accordance with the Condon approximation, and |F,/), the unnormalized vibrational wave 
function of the excited molecule. Due to (1) this wave function consists of two terms 

\Vj)=\v,j) + a\Vj), (11) 

, e functions \ v. ,/) correspond to the vibrational wave packets being excited respectively by 

each of the subpulses in (1). Generally, the wave functions | V ,/) and \ v ,/) are not orthogonal and 
their overlap depends on the delay time T between the subpulses. This opens a possibility to 
observe the delay time dependent interference effects in the physical processes which are 
determined by the population of the excited electronic state of the molecule. For the sake of 
explicitness we shall assume that the measuring quantity is the quantum fluorescence yield of the 
excited molecule. In the case of the transition in the pure electronic system (without a vibrational 

degree of freedom) functions \ v, ,/) are reduced to numerical amplitudes 


V , = 


n 


exp 


8 2 _ iST 
+ 


2 u 2 ' 2 

where 8 = s e jh- ejh- Q. The dependence of the excited state population on the delay time 
between the subpulses in this case is of a trigonometric character: 


( 12 ) 
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(13) 


In the remainder part of the paper we shall consider how the accounting of the molecular 
vibrations changes Eq. (12). The general expressions for functions \v ± ,t) are given by 


K>') = ““T~ J dt ■ ex P ” 1 T A)~ ,Q/ ' + ~ ^ + J £e ‘ 
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Using (7) and (9) it is convenient to transform the wave functions (14) into the form of a 
distributed coherent state [13] 
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dt x exp 
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Time dependence of the Heisenberg displacement operator D(l ) in (15) is determined by the 
Hamiltonian (6). To evaluate the time integral in (15) now it is possible to use the well known 
Fock-space expansion of the coherent state £>|0) [23]. As a result we have: 
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(16) 


|/t,/) = exp(-/«ft>/)|«), 

where |/t) are the Fock-space eigenfunctions of the harmonic oscillator. Comparing Eqs. (12) and 
(16) it is easy to see that the coefficients at 1 0, /} in (16) coincide with the amplitudes V , . Exact Eq. 
(16) is convenient for the numerical calculations. Another way to evaluate the integral (15) is based 

of the coherent state £>(' ~ ^ 


on the replacement 


it, + y 2 -/)|o> = 
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by 
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-ge ' 4 * r ’>0 -icot^fj, which is possible when u»co. In [13,15] it was shown that the integration 


in (15) in this case can be performed in the operator form, and leads to the result: 


V ± ,t) * W ' 2 exp 


(_ .ST 


&{' ± A) s ( i ± T A )\ °>- 


(17) 


In (17) the normalizing factor W is the absorption spectrum of the electronic- vibrational 
transition: 


1 u 

W = n 2 ——exp 
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(18) 


D and S are the unitary displacement and squeezing operators: 


D(0 = expi-g 


1 


co A 

~bJ 


[/> (/) -/>(,)] , 
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S(t) = expj —In — I [/>(0 2 - (/) 2 ] 


V u 


b(t ) = exp(-icot)b, 


(19) 


The action on the vacuum state of the operator DS gives an ideal coherent squeezed state in which, 
generally, the mean position and momentum of the oscillator have nonzero value, the uncertainties 
of the position and the momentum are not equal to each other and their product has the minimal 

possible value. In the a-plane of the coherent states each of the vectors | V, ,/^ can be considered as 

an ellipse uniformly moving along the circle with the radius g = an< ^ keeping the 

orientation of the axes in the rotating frame. At moments t , for which the arguments of the 
operators in (17) are equal to the integer number of the vibrational periods, the values of the ellipse 
main axes are equal to the non-dimensional uncertainties of the position and the momentum with 
the small axis corresponding to the uncertainty of the position: 

A? = F i i , =t * J2g 2 ty 2 +h 2 

Ap = — —y = . (20) 


V2-y/2g 2 iy 2 +M 2 


■Jlu 


If the coupling constant is large enough (g>>l), then the wave functions are located at a 

sufficient distance from each other, and are approximately orthogonal except for the case when the 
time delay between the subpulses appears to be equal to an integer number of the periods. From 
geometrical considerations it is clear that when the delay time between the subpulses gains an 

integer number of the vibrational periods, the overlap of the vibrational wave functions | F + ,f) and 

| V ,t) can be observed in the a-plane as the approaching of two ellipses along the direction of their 
big axes. So it is likely to expect the following qualitative effect of the molecular vibrations on the 
interference picture given by Eq. (13): (i) the dependence of the population of an excited state on 
the delay time between the subpulses should be more pointed in comparison to the trigonometric 
dependence in Eq. (13), (ii) while the delay time approaches the integer number of the periods the 
interference picture should be sensitive to the quantum state of the molecular vibrations, and the 
squeezing effect can be seen due to the dependence of the interference picture on the subpulse 
duration. 
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Let us put r = + r in Eqs. (15), (16), where r , is a small deviation from the integer number 

of the vibrational periods, and evaluate the population of the excited state. In accordance with 
Eqs. (15) and (16) one finds 


(V\ V) = u] 
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drjexu iSrj - ~T r l‘ 
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( 21 ) 
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(22) 


where M(^) = ^o|z> + (0)Z>(^ )|o^, and 

/ , v i- v^g 2 " r {s + nco) 1 
(V\v) = 2n 2 SVexp -g 2 - 

While deriving Eq. (21) we have performed the integration over the summary time |(f, +r, ). The 

remaining integration in (21) is carried out over the difference time t x -t x =tj. The vacuum average 
of two displacement operators in (21) frequently arises in the theory of Franck-Condon transitions 
and can be evaluated, for example, using the coherent states method [23]. 

M(£)=ex p[g 2 l)]* ( 23 ) 

It is worth noting that Eq. (22) can be received from (21) by expanding in (23) the exponent in the 
powers of g 2 exp(/<y£). Due to the condition u»co integral in (21) can be evaluated 
asymptotically (by the steepest descent method [24]). To do that one can expand the index of the 
exponent in the powers of % up to the second order terms: 

M{£>) = exp(/g 2 <u£ -yg 2 <y 2 ^ 2 ). ( 24 ) 
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(25) 


equation UcSCIlUCS me linciititiivt uunvwi uiv , ^ 

molecule. Depending on the experimental conditions, i.e. the values of cp , S, and m, this 
interference can lead to a sharp peak (constructive interference), or to a down-fall (destructive 
interference) in the fluorescence dependence on the delay time between the subpulses. The new and 
the most constitutive feature of Eq. (25) is the dependence of the interference peak width y 0 on the 
reciprocal pulse duration u. Within the range of pulse duration (a 2 «u 2 «2g 2 a> 2 , where, in 
accordance with Eq. (20), the squeezing effect is the most considerable, the width of the 
interference peak approximately equals to the subpulse duration y 0 . This gives a practical 
possibility to use this effect for the duration measurement of a femtosecond laser pulse. The fact 
that the dependence of the interference peak width on the pulse duration is indeed the consequence 
of the dynamical squeezing effect becomes clear if we note that the second order expansion in Eq. 

(24) is equivalent to the first order expansion over /, ± x /i * n Eq. (15). In the limit of extremely 
short pulses w 2 » 2 g 2 co 2 Eq. (25) coincides with that previously received in [20]. In that limit the 
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interference peak width does not depend on the pulse duration: y 0 %(V2g<y) ' .It should be noted, 

that the limiting interference peak width is achieved for those laser pulses, during the action of 
which the initial vibration state of a molecule in the process of electronic excitation has no time to 
change. The effect of the dynamical squeezing in this limit is absent. The picture of the interference 
peak as the function of the variables r and u , numerically evaluated with the help of Eq. (22), 
can be seen in Fig. 1 . 


The case of nonzero temperature. In the case of nonzero temperature we shall assume that the 
initial vibrational state of the molecule is described by the equilibrium density matrix. The 
population of the excited electronic state in this case can be evaluated in the second order 
perturbation theory for the density matrix. An expression for the nonequilibrium density matrix of 
the phonons, arising due to the excitation of a Franck-Condon transition by a single laser pulse 
was obtained in [11, 14]. In the present paper we need to calculate the trace of an analogous 
density matrix, appearing due to the action of the field (1). It seems to be evident, that the 
expression for this trace will have the form of Eq. (21) with the functions M{4) replaced by 

(M^)) = ^[ A Z) + (0)^)] = exp{^[(/t + iy^ +he -2n-l]j, (26) 

where p 0 is the equilibrium density matrix of phonons, and n, the thermal equilibrium phonon 
number. The Fourier expansion of (a/(^)) is given by 


(M^)) = Z I' ^-J A,(V 4n{n + 1) ) exp[-g 2 (2/i + 1) + in 6)%] , (27) 

h--co v ” y 

where is the modified Bessel function. A generalization of equation (22) for the case of a 
nonzero temperature has the form 

(V\V) = 2 tt' 2 ^(W"(« + 0)exp -g 2 (2n + l)- 

n=-<x\ " y 


(S + n co ) 3 


r s 

l + |a| + 2|ar | cos|^ + 2n — + (<5 + nry)r 


(28) 


To evaluate the integral in Eqs. (21), in analogue to the case of the zero temperature, we shall 
expand the index of the exponent in Eq.(26) in the powers of £ up to the second order terms. In 
this way we find that the effect of finite temperature in Eqs. (18) and (25) is manifested in the 
replacing coefficients B 0 and y 0 by 


B= yju 1 + 2g 1 <v 1 (2n + 1) , 


^u 1 +2g 2 <y 2 (2/J + 1) 


UgCOyj 2( 2/1 + 1) ^ ) 

From (29) one can see that the temperature growth leads to the extension of the area of linear 


dependence of the interference peak width on the laser pulse duration. According to. Eq. (29), 

w <> X =W 0) 2 (2h + \)Y 2 can be considered as the limiting value of pulse duration for the range of 
this linear dependence. It is interesting to estimate the value of u Q 1 for a real physical system. To do 


this we shall use the data for I 2 molecule from [22]. Taking g 2 = 6.4, = 300/s and T = 300°/^ , 

we obtain w 0 ‘ * 7 fs. The natural upper limit for is the inverse vibrational frequency co In the 
experiment [22] the laser pulses of about 50 fs duration have been used. In accordance with the 
given estimate the pulses of such duration belong to the area where y « w '. 
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In conclusion note that the measurements of the interference peak width aimed at determining 
the duration of a femtosecond laser pulse represents only one of the possibilities to use the 
intramolecular interference phenomenon of the vibrational wave packets. Another interesting and 
important possibility is the observation of the wave packet distortions due to the intramolecular 
propagation. This problem we are planning to consider in our forthcoming paper. 
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